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One of the tasks of urban and hazard planning is to mitigate the damages and minimize 
the costs of the recovery process after catastrophic events. The rapidity and the effi- 
ciency of the recovery process are referred to as resihence. A mathematical definition of 
resilience of an urban community has not yet been identified. In this paper we propose 
and test a methodology for the assessment of urban resilience to a catastrophic event. 
The idea is to merge the concepts of engineering resilience and ecosystem resilience. As 
first step we model a urban community by means of different hybrid complex networks, 
composed by human elements, i.e. the citizens, and physical elements, i.e. urban life- 
lines and infrastructures. As second step, we define a class of efficiency indexes on these 
hybrid networks. Then, we quantify the resilience of a urban system by computing these 
efficiency indexes before a simulated catastrophic event, after the event has occurred and 
during the subsequent recovery process. As a case study, we test this idea in case of 
simulated earthquakes in the city of Acerra, Italy. 
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I. INTRODUCTION 

Resilience of urban systems is a key issue for contempo- 
rary society. Natural and human- induced risks threaten 
contemporary cities more than ever, because of the in- 
creasing urban population. The intensity of natural ex- 
treme events is in many cases also exacerbating, because 
of climate change phenomena. Urban management needs 
to cope with this issue and to implement risk mitigation 
and resilience strengthening policies. The assessment of 
city resilience against natural and human-induced haz- 
ards represents an urgent challenge for scientific commu- 
nity. This has been already highlighted by D alziell and| 
(2004|), which pointed out the need to intro- 
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duce appropriate metrics to quantify resilience. 

However, a unique definition for urban resilience 
against extreme events is not available. So far, two 
philosophies have inspired the recent approaches assess- 
ment the resilience of urban systems against disasters: 
the vision and the ecosystem vision. According to the for- 
mer, city resilience depends on the capability of all the 
physical components, to absorb an external shock and 



quickly return to the initial state (Rourke, 2007 Reed 
eFldl [20071 IBruneau et all [2003) |Pimni| p^Ml [OprT 



covic and Tzeng, 2002). The engineering approach can 



be straightforwardly used to quantify resilience of single 
urban components/infrastructures, by quantifying their 
time- varying performance in the aftermath of the shock. 



In the ecosystem approach, instead, resilience is defined 
as the capability of the whole urban system to recover 
with the full functionalities that existed before the shock. 



even without returning to the initial state ( Rolling, 1973 



1986, 2001 Kovalenko and Sornette, 2012). This implies 



that the system may evolve in something slightly dif- 
ferent, keeping its identity in a broader sense, as well 



described by the theory of complex systems (Bar- Yam 



2003). 



Complex systems are often defined as systems consist- 
ing of many interacting units with the ability to generate 
a non-trivial, collective behaviour through the combina- 
tion of simple mechanisms acting at a local scale. In such 
systems, the evolution of the whole is not readily pre- 
dictable by observing the behaviour of its components. 
From this point of view, cities and urban systems are 



undoubtedly complex systems (Batty 2005). They have 



evolved over the centuries from the bottom up before 
reaching the actual configuration, under the concurrent 
action of a spontaneous unsupervised process, which has 
continuously shaped the urban structure at a microscopic 
level, combined with more recent planning efforts, which 
have tried to force more efficient structures on top of 
a pre-existing urban backbone. In fact, a city can be 
thought as a system composed of at least three layers. 
The first one is the technological layer, which consists of 
the buildings, the services and the infrastructures that 
physically link them together. The second one is a ge- 
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ographical layer, which takes into account the fact that 
cities are embedded in a geometrical space and their de- 
velopment is affected by the surrounding environment. 
Finally, there is an human component, i.e. the social 
network of citizens that reside, move and work within 
the physical frame. 

Apparently, the ecosystem approach is more suitable 
for urban systems. Nevertheless, the engineering ap- 
proach is also extremely meaningful, since the engineer- 
ing resilience holds in the capability to reach a new state 
that is not "worst" than the pre-shock one. Hence, 
metrics to quantify the pre-shock and post-shock perfor- 
mance of a city and to assess whether the pre-shock per- 
formances are recovered after the event are fundamental 
tools to analyse the resilience of a city in a quantitative 
way. In this article we propose a series of measures to 
estimate the pre-shock and post-shock performance of a 
city, as required by the engineering approach, which take 
into account the possibility that the urban system evolves 
in a different configuration, as allowed by the ecosystem 
approach. We have also used these measures to test and 
compare three different strategies of reconstruction. 



II. A METHODOLOGY TO QUANTIFY CITY 
RESILIENCE 

The methodology presented here aims at identifying a 
unique approach to quantify, from an engineering point 
of view, the resilience of a city. However, inspired by 
the ecosystem vision, it moves from the point that urban 
systems can have a resilient response to a shock, also at- 
taining a new configuration of the physical systems, that 
is a new dynamic equilibrium stage. In other words, we 
believe that a city can be considered resilient although 
some of the infrastructures do not recover their full func- 
tionality after the shock and new infrastructures are built 
elsewhere to recover the performance of the whole urban 
system. Given this principle, a first issue is to identify 
a proxy of the global urban system performance, to be 
computed before and after the shock and during the re- 
covery phase in order to quantify the capability of a city 
to recover its performance before the shock occurred. We 
move from the observation that the street network is one 
of the most important aspects of a city, since its struc- 
ture is intimately connected with the reachability of ser- 
vices and facilities and therefore with the overall quality 
of life perceived by the citizens. Consequently, we con- 
sider the efficiency of the street network of a city, appro- 
priately weighted to take into account local population 
density and the physical location of services and facili- 
ties, as a meaningful proxy of the efficiency of the city 
as a whole. We make use of some basic concepts and 
metrics of complex network theory in order to quantify 
this efficiency and to compare the efficiency of different 
configurations. Complex network theory has proven to 



be a robust theoretical framework to study the struc- 
ture of networked systems and has been largely employed 
in the last decade for the characterization of a variety 
of phenomena occurring in systems composed by inter- 
connected units, including many biological, technologi- 



cal and social networks (Strogatz 2001 Newman 2003 



Boccaletti et al., 2006). Recently, the complex networks 
approach has also been successfully employed to quan- 
tify and model the structural aspects of spatial networks 



([Barthelemy 


2011) and of street networks 


(jCrucitti et al. , 


2006 Strano et al. , 


2012). 



A. Networks of urban street patterns 

The functioning of the infrastructures and services of 
a urban agglomerate heavily relies on the existence of 
an underlying street network, and their efficiency usu- 
ally depends on the structural properties of the urban 
street network. A huge amount of literature in urban 
morphology indicates that roads are a fundamental driver 
in urban evolution ( Marshall) [2006 Southworth and Ben^ 
Joseph, 2003) and complex network theory has provided 



significant contributions to the quantitative characterisa- 
tion of urban street patterns, as reviewed in |Barthelemy| 
(2011). More recent studies have also shown that road 



networks are important for the dynamical processes oc- 
curring on them and for the evolution of urban areas 



in general (Batty, 2005 Strano et al. , 2012 Bettencourt 
et al.[ |2007[ [Balcan et al.[ |2009[ ). Reliability of infras- 



tructure networks prone to natural hazards is largely dis- 
cussed in available literat ure (jPinto et al. 2006). Con- 



nectivity based methods (Li and He 2002 Duehas Os- 



orio and Rojo, 2011) and flow based methods (Franchin 
and Cavalieri, 2013) have been proposed and discussed 
to assess the reliability of lifelines in damaged config- 
uration. Generally, networks can conveniently be de- 
scribed by means of graphs consisting of a set of points 
A/", called nodes or vertices, and by a set V of edges 
connecting pairs of points. A graph with N = \JV\ 
nodes (A/" = {ni, n2, na, . . . nAr}) and K = \V\ edges 
(V = {vi^V2^ vs, . . .Vk}) can be represented by giving its 
adjacency matrix, i.e. the N x N matrix whose entry aij 
is equal to 1 if there is an edge connecting node i and 
node j, while a^j = otherwise. It is also possible to 
assign a weight or a length lij to each edge linking nodes 
i and j in a graph, thus defining a weighted adjacency 
matrix {hj}. 

Spatial networks are a special class of complex net- 
works whose nodes are embedded in a space associated 
with a metric. Typical examples of spatial networks 



include electric power grids ( [Kinney et al. , 2005) and 
transportation systems including rivers, trade routes and 



street networks (Crucitti et al. 2006 Strano et al. , 2012 



Pitts 1965). In the case of street networks, each cross- 



ing is represented by a node while edges represent street 



3 



segments, so that two nodes are connected by an edge if 
the corresponding crossings are adjacent to the same seg- 
ment of road. In the fohowing we denote by G(A/', V) the 
graph representing the street pattern, where M is set of 
street crossings and V is the set of street segments. Street 
networks are naturahy embedded in a two-dimensional 
Euchdean space, whose metric is the usual Euclidean 
distance, so that the lengths Uj of the edges satisfy the 



triangular equality (Barthelemy 2011). 



In transportation and communication networks it is 
usually important to know how to move or send an in- 
formation from a node i to another node j. An alternate 
sequence of nodes and edges that starts from i and ends 
in j is called a walk from i to j. If there exists a walk 
between node i and node j, we say that i and j are con- 
nected. A maximal set of nodes which are mutually con- 
nected to each other is called a component of the graph. 
If not all the pairs of nodes in the graph are connected, 
then the graph is composed by more than one component. 
Each walk is associated to a cost, that is the sum of the 
lengths of the edges involved in the walk. If each node of 
the walk is traversed only once, then the walk is called 
a path. The path from i to j havbning minimal length is 

If 

is set to 

oo, and the two nodes are said to be disconnected. 

A measure of the typical separation between nodes in 
the graph is the characteristic path length L, that is the 
mean value of the length of the shortest paths between 
all the possible pairs of nodes: 



called shortest path and itsa length is denoted by dij. 
two nodes are not linked by any walk, then d^j 



N(N -1) 



(1) 



In general, the lower the characteristic path length, the 
better the communication between any pair of nodes cho- 
sen at random. Using the characteristic path length to 
measure the resilience of a city is possible but indeed not 
convenient, since L becomes infinite as soon as there ex- 
ists at least one pair of disconnected nodes. However, 
the result of a shock event on a city, like an earthquake 
or a fiood, is often a disconnected street network, which 
always has an infinite characteristic path length indepen- 
dently of the actual number of pairs of sites that remain 
disconnected after the event. The efficiency of informa- 



tion exchange over a network, proposed in reference ( La- 
[tora and Marchiori , 2001), is a measure which allows to 
overcome the subtleties due to infinite characteristic path 
lengths and can be therefore used to quantify the aver- 
age reachability of the nodes of a network after a shock 
event. The efficiency e^j of the communication between 
nodes i and j in a generic graph is deffned as the inverse 
of the length of the shortest path connecting i to j, i.e. 
Cij = 1/dij. The efficiency is minimal when i and j are 
disconnected, i.e. when dij = oo. In the case of spatial 
graphs, it is possible and convenient to normalize the 



efficiency of a pair of nodes by the Euclidean distance 
between them, defining e^j = dfj'^^/dij, so that the re- 
sulting normalized efficiency is maximal and equal to 1 if 
and only if the shortest path between i and j runs exactly 
along the direction of the geodesic which connects them. 
The global efficiency of a spatial network is defined as 
the average of the normalized pairwise efficiency over all 



possible pairs of nodes (Vragovic et al. , 2005|: 



E 



N{N- 



^eucl 
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E 



(2) 



Notice that the global efficiency is normalized in [0,1], 
since in general the distance dij between node i and node 
j, which is measured as the total length to be traversed 
in order to get from i to j using a sequence of street 
segments, is larger than the Euclidean distance between 
i and j, so that each term in the summation is < 1. 
Consequently, it is possible to compare in a consistent 
way the efficiencies of two distinct graphs G and G' , even 
if they have a different number of nodes and edges. 



B. Hybrid social-physical networks, HSPNs 

The main idea behind the proposed methodology is to 
define Hybrid Social-Physical Networks (HSPN), and to 
assess their efficiency, before, during and after an extreme 
event. Basically, HSPN are physical spatial networks, 
composed by infrastructures and facilities, enriched and 
enlarged with further nodes, representing the citizens, 
and with further links, representing the interactions be- 
tween citizens and physical components. One HSPN can 
be defined for each of the physical networks interacting 
with citizens. More specifically, one HSPN can be identi- 
fied for each service provided by urban physical systems 
to citizens. For example, a HSPN can be defined for res- 
idences, namely the residential HSPN. In this case two 
sets of nodes and two set of links are identified: 

• The set of intersection nodes, which we called M. 

• The set of building nodes, whose elements are the 
centroids of residential buildings (hereafter referred 
as B). 

• The set of street segments, which we called V. 

• The set of door links, whose elements represents 
the logical links connecting a building node to an 
intersection node, as showed in figure [l] This set is 
donoted by Vb 

The street network augmented with the set of build- 
ings and door links is denoted by G23(A/' U S, V U V23). 
This network is enlarged with the social components to 
build the HSPN: each residential building is linked to 
further ending nodes, representing the citizens living in 
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FIG. 1: Network representation of a city. Buildings 
nested in the same point have nuh distance. The nodes 
representing these buildings belong to the set X. 



of that service by the citizens, i.e. as a measure of the 
quality of the service provided. Thus, it is possible to 
quantify and compare the performance of the HSPN in 
distinct configurations of the physical networks, even if 
some infrastructures are not present in one of the configu- 
rations. By performing this analysis before and after the 
reconstruction which follows an hazardous event, which 
damage some of the pre-existent physical components of 
the urban system, we can quantify how the quality of 
public services provided to citizens has been affected, al- 
though the physical systems are rebuilt and rearranged 
in a new configuration. The idea is that this could rep- 
resent an effective methodology to quantitatively assess 
the resilience (as suggested by the engineering resilience 
approach), but with a systemic perspective and focusing 
on the social components of cities (as indicated by the 
ecosystem resilience approach). 



each building. The new links represent the interaction 
between citizens and buildings, that is the interaction 
between the social components and the physical compo- 
nents. This is clarified in Figure [T] The obtained net- 
work can be analysed in terms of different network per- 
formance metrics, defined in section [II. C , to appreciate 
the capability of the urban system to provide services to 
citizens, that is, in this case, to allow citizens to live and 
interact with each other. Both the social and the physi- 
cal part of the HSPN can be even more complex. For in- 
stance, if we consider the street network, toghether with 
residential and commercial buildings, we can estimate 
the capability of citizens to reach goods supplies. The re- 
sulting graph is denoted by G^(A/' U S U 0, V U U V^), 
where Q is the set of nodes representing commercial 
buildings and Vq is the set of logical links connecting this 
commercial buildings with the adjacent s street segments. 
In graph Qq the citizen nodes are attached to their cor- 
responding building nodes, while goods nodes are con- 
nected to the commercial buildings where they are avail- 
able. In the following we refer to this augmented network 
as to the goods HSPN. Similarly, a service HSPN can be 
built to assess the capability of citizens to be linked to 
schools, hospitals and other public infrastructures and to 
receive public services. In this case the underlying phys- 
ical network is denoted by G5(A/' U S U 5, V U U V5), 
where S is the set of nodes representing public build- 
ings and Vs is the set of logical links connecting public 
buildings to intersection nodes. 

Other lifelines can also be represented by means of ap- 
propriate HSPNs. For instance, the electric grid network 
and the water supply /sewerage network can be intro- 
duced, respectively, to investigate the capability of the 
urban system to provide citizens with electricity and wa- 
ter. 

The potential of this approach lays in the fact that the 
global performance of the HSPN associated to a given 
service can be considered as a proxy of the accessibility 



C. Measuring the performance of HSPNs 

We define here some efficiency metrics for HSPNs 
which are inspired by the normalized global network effi- 
ciency given in equation ([2|. The mere definition ([2| is 
inadequate to estimate the efficiency of the city, because 
it does not take into account the number of people liv- 
ing in each building. We can then devise a measure of 
efficiency of communications between people inside build- 
ings by summing over all couples of inhabitants and using 
df^^^ /dij = 1 for couples of inhabitants living at distance 
zero, i.e. assuming the maximum efficiency in their com- 
munications. Finally it is possible to calculate an expres- 
sion of efficiency of communications of citizens, similarly 
to equation ([2|, as: 



Htot{Ht 



tot 



■ E 



^eucl 



(3) 



where z, j are the indexes of nodes representing buildings, 
Hfot is the total number of inhabitants of the city. Hi is 
the number of people living in building z, S is the set of 
nodes representing buildings and dij is the length of the 
shortest path between i and j evaluated on graph G^- 
Notice that in the summation over j, j ^ i we assume 
that dfj"^^ /dij = 1 for couples of buildings at zero dis- 
tance (linked to the same street node). The term {Hi — 1) 
inside the round brackets, times i^^, takes into account 
the couples of inhabitants that live in the same building 
and whose efficiency is en = 1. We choose this measure 
to estimate the global efficiency of citizen-citizen commu- 
nication in a city. 

It is also possible to rewrite equation ([3| in the follow- 
ing form, in order to avoid to impose d^^^^ / dij = 1 for 
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couples of buildings at distance zero: 



Htot{Ht 



tot 



E 



H. 



^euc 



(4) 

where Htot is the total number of inhabitants of the city, 
Hi is the number of people living in building i and hi 
is the number of inhabitants that live in the set X of 
buildings with zero distance to building i; the summation 
over j is now limited to couples of nodes with non zero 
distance. 

It is possible to write similar expressions also for goods 
HSPN by summing the outer summation over the set Q 
of the buildings that contain goods, for instance shops, 
and normalizing by the total number of goods in the city 
Gtot^ thus leading to equation: 



E. 



Crtot^tot ^^^^^^ 
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(5) 



dij 



where Gi is an estimate of the amount of goods in the 
shop i G ^, dij is the length of the shortest path between 
i and j evaluated on the graph Gg . 

Similarly, an estimate of the performance of citizens- 
schools HSPN is given by equation 
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S tot Hi 
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^eucl 
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where Si represents a measure of the importance of the 
school i (e.g. the floor area of the building or the number 
of students), S is the set of nodes representing school 
buildings and Stot = Xli '^i- ^^^^ equation the length 
dij of a shortest path is evaluated on the graph Gs- 

The main idea behind this approach is to use these 
measures as a proxy of the quality of the relationships be- 
tween people and infrastructures; that is, in other words, 
through these measures one can appreciate how services 
of physical systems are available. Notice that, in the 
HSPNs efficiency quantification the importance of the fi- 
nal nodes of the physical network (i.e. the buildings) 
needs to be "weighted". To do this, the floor area was 
used; in fact, it can be interpreted as a proxy of the num- 
ber of citizens living in each building, for the residential 
HSPN, of the amount of goods in each shop, for the goods 
HSPN, and of the importance of each school building, for 
the for citizens-schools HSPN. 




FIG. 2: Acerra location 




FIG. 3: Acerra view 



III. A CASE STUDY: THE CITY OF ACERRA 

Acerra is a town in the Province of Naples, in Italy, 
about 20km north-east of Naples (Figure |2|; it is one of 
the most ancient cities of the region, probably founded 
before 400 BC. The urban configuration, commonly expe- 
rienced in medium cities in Italy, is composed by a dense 
historical centre, with mainly ancient masonry buildings, 
surrounded with more recent urban expansion areas, with 
mainly reinforced concrete buildings; the whole built-up 
area is surrounded by a countryside area, where a large 
industrial area is also located (Figure |3|. The whole ter- 
ritory is extended over a surface of about 54km^ and its 
population is estimated about 55000 inhabitants. 

Acerra is prone both to seismic risk and flood risk. In 
particular, the seismic risk is generated by the seismo- 
genetic areas located in the Appennines zones, at nearly 
one hundred kilometers; Figure [4] reports the seismic haz- 
ard in terms of annual rate of occurrence of events with 
Pick Ground Acceleration (PGA) larger than a certain 
value. Flood risk is due to the Regi Lagni river, bor- 
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FIG. 4: Seismic hazard: annual rate of exceeding PGA 
vs PGA 



dering the built-up area. Furthermore, an industrial risk 
can be also identified, because of the infrastructures and 
factories located in the industrial area. 



A. Quantifying city resilience 

The objective of this section is to practically deploy the 
proposed methodology, focusing on the seismic risk. The 
final aim is to quantify the aptitude of the city of Acerra 
to exhibit a resilient response, in case of earthquakes. To 
do this, the following steps have been implemented. 

1. Data related to the buildings and the street pat- 
terns are collected by means of a GIS support. In 
particular, the buildings are divided into structural 
typologies and the total fioor area is measured for 
each building. Data related to each street link is 
also collected. The network composed by street 
patterns and building is depicted in figure [5] 

2. A fragility model is assigned to each building, given 
the structural typology; in other words, for each 
PGA value, the probability that the building over- 
come a specific damage state has been defined. In 
details, the "onset of damage" limit state has been 
considered, as the damage state making the build- 
ing unfit for use, whatever is the use of the build- 
ing. Furthermore, a damaged building can also in- 
terrupt the transit along the streets ahead. Thus, 
the probability that the building, once beyond the 
"onset of damage" , could make the street unfit for 
passage, was also introduced. It was made propor- 
tional to the ratio of the height of the building over 
the width of the street. 

3. In each residential building, the number of residents 
is estimated. A similar estimation on the impor- 
tance of the buildings is made for other non resi- 



dential buildings. In example, for the commercial 
buildings and the school buildings, an estimation of 
the amount of goods and the number of students, 
respectively, is needed. Hence, different HSPN are 
built, joining physical networks, made by buildings 
and street patterns, and social networks made by 
residents, goods, students, etc. 

4. Given the fragility model, different earthquake sce- 
narios are generated, characterized by a certain 
value of PGA. In particular, in each scenario and 
for each building, the occurrence of a damage state 
beyond the "onset of damage" limit state is veri- 
fied, by means of a Monte Carlo simulation, given 
the PGA value and the fragility model. Similarly, 
for the damaged buildings, the occurrence that the 
streets ahead are interrupted is verified by means 
of a Monte Carlo simulation, given, in this case, the 
probability that the street is made unpracticable. 
Then different recovery strategies are considered, 
each of them consisting of different urban configu- 
rations, where part of the damaged buildings and 
of the interrupted street patterns are progressively 
recovered, up to the reallocation of all the displaced 
citizens. 

5. Finally, the generated HSPN are characterized by 
means of efficiency measures that are computed be- 
fore and after the earthquake scenarios and after 
each reconstruction step of the considered recovery 
strategies. 



B. Numerical analysis and results 

We have implemented the aforementioned procedure 
for the city of Acerra in different earthquake scenarios, 
ranging between 0.05^ and 1.0^. We assume that the 
fragility models providing the probability of exceeding 
the "onset of damage" limit are characterized by a log- 
normal distribution function for both masonry and rein- 
forced concrete buildings: 



P5 (x;/i,cr) 




(7) 



where erfc is the complementary error function and the 
variable x represents the value of PGA; the parame- 
ters /i, cr have been considered equal to —1.03,0.35 for 
masonry buildings and 0.91,0.29 for reinforced concrete 



buildings, according to reference (Ahmad et al. , 2011). 
Furthermore, it was supposed that a road in front of a 
damaged building is made impracticable with probability 



p ^) _ J 1 iff > 

' 1 J otherwise, ' 



(8) 



FIG. 5: Buildings (black points) and street patterns (red lines) network for the city of Acerra 



where h is the height of the building and / is the width 
of the road. To estimate the number of residents in each 
building, one inhabitant per 30 square meters was as- 
sumed. This value was obtained by dividing the total 
number of inhabitants, i.e. 55000, by the total floor area 
of all the residential buildings. HSPN were built also 
for goods and students. In these cases, the floor area of 
commercial buildings and schools was used, respectively, 
as a proxy of the importance of each building. A set of 
measures have been carried out in the damaged config- 
uration, after the considered earthquake scenarios and 
the values have been compared with those obtained in 
the undamaged configuration. In particular, three sets 
following measures have been computed: 

1. The first set is composed by the efficiencies Ecc^ 
Ecs and Ecs- The obtained values are normalized 
by the corresponding values in the undamaged con- 
figuration and are plotted vs. the PGA in figure [7| 

2. A second set is composed by the following mea- 
sures, whose values are reported in figure [6] vs. the 
PGA: 

• The number of undamaged buildings; 

• The number of not displaced citizens; 

3. A third set of measures is performed on the net- 
work of street pattern G defined in section II. A 
Let us recall that a street in front of a damaged 
building is feasible to be blocked. We simulate this 
process by removing the intersection nodes linked 
to a damaged building with probability P^. The re- 
moval of such nodes implies also the removal of the 
incident links. The deletion of links could result in 
the fragmentation of the network in several com- 
ponents. This possibility could be disastrous, since 
it some parts of the city will remain isolated from 



TABLE I: Absolute values in the undamaged 
configuration 



Ecc 0.746 

Ecg 0.252 

Ecs 0.313 

undamaged buildings 3493 

not displaced citizens 55000 

Number of connected components 1 

Characteristic path length (m) L(0) = 599, 2 
Nodes belonging to the largest component S{0) = 4638 



the others. Then, it is possible to study how ex- 
treme events damage the street pattern measuring 
the following quantities: 

• The number of connected components; 

• The number S of intersection nodes belonging 
to the largest connected component; 

• The characteristic path length L of the largest 
connected component, as defined in eq. Q; 

In the figures [7| [6] and |8] the values for the undam- 
aged configuration are reported in correspondence of zero 
PGA. The absolute values are also reported in table [Tj 

It can be observed from figures [6] that the number of 
undamaged building and the number of not displaced 
people have an exponential decay with the PGA with 
exponent b —4. We observe a faster exponential decay, 
with exponent b —9, for Ecc, Ecg and Ecs for values 
of PGA in the range between 0.25 and 0.6, as showed 
in figure [7[ These results confirm that the efficiency of 
HSPNs can be reliably used as a proxies of the quality of 
an urban system. 
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FIG. 6: Undamaged buildings and not displaced 
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FIG. 7: Efficiencies Ecc^ Ecg and Ecs vs PGA over the 
pre- shock values. 



Furthermore, from figure [8] we note that behaviour of 
the street network is consistent with a percolation tran- 
sition, indicating the existence of a critical range around 
PGA = 0.3 — 0.4 beyond which the street network is bro- 
ken into many parts and loses the giant connected com- 



ponent (Stauffer and Aharony , 1994; Dorogovtsev et al. 



2001 ). This is due to the adopted fragility functions that, 
in correspondence of this range of values, present a sig- 
nificant increase of the probability of collapse. 

Reconstruction is analysed only for 3 earthquake sce- 
narios, namely 0.2, 0.25 and 0.3 PGA. The configuration 
of the damaged physical network after each scenarios is 
reported in Figure [9j For each scenario, 4 reconstruction 
strategies are considered. 

The first one, hereafter referred to as status quo down- 
up, consists in striving to return the city to the pre- 
event configuration. The process of reconstruction is dis- 
cretized into n steps. In each step a fraction 1/n of the 
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FIG. 8: Street network fragmentation under earthquake 
event. The number of connected components (red line), 
the normalized characteristic path length of the largest 

component L{PGA) / L{fi) (green line) and the 
normalized size of the largest component S{PGA)/S{{)) 

(red line), as function of increasing intensity of 
earthquake, measured in unit of PGA. 



displaced citizens (those living in the damaged buildings) 
is allocated, assuming to rehabilitate the buildings from 
the smallest to the largest, i.e. from the cheapest to 
the most expensive. Blocked roads are recovered when 
the buildings that caused their interruption are recon- 
structed. In the second strategy, hereafter referred to 
as status quo up-down, as in the previous one, the city 
returns to the undamaged configuration, but the build- 
ings are restored from the largest to the smallest. The 
third strategy, hereafter referred to as new sites down-up, 
consists of reallocating part of the displaced citizens in 
new residential sites. Also this process is discretized into 
steps. In the first step the new buildings are realized in 4 



empty areas (highlighted in figure 10), to reallocate the 
20% of the displaced citizens. Then the existing build- 
ings are rehabilitated as in the first strategy, from the 
smallest to the largest ones, until all citizens have been 
reallocated. This strategy is often adopted in real cases 
to recover as soon as possible urban functionalities and 
has been experienced in several disaster recovery strate- 
gies, e.g. in the aftermath of the L'Aquila 2009 earth- 



quake, in Italy, where the CASE project (Cosenza and 
Manfredi, 2010[ ) was implemented. In this strategy some 
of the existing buildings are not recovered, since a part of 
the population is reallocated elsewhere. To re-establish 
the original urban street pattern, it has been assumed 
that the interrupted roads, that would have not been 
recovered (since interrupted by buildings that are not re- 
covered), are re-established at the last step. The fourth 
strategy, hereafter referred to as new sites up- down, is 



(a) 



(b) 




(c) 

FIG. 9: Configuration of t lie damaged piiysical network after tiie 0.2PGA(a), 0.25PGA(b) and 0.3PGA{c) 

earttiquake scenarios. 



similar to tiie previous one, but, after reallocating the 
20% of the citizens in the new buildings, the existing 
buildings are restored form the largest to the smallest 
ones. The last two strategies, namely status quo inwards 
and status quo outwards consist of rebuilding the city as 
in its undamaged configuration, moving from the sub- 
urbs to the center of the city and viceversa, respectively. 
After each step of the reconstruction strategies the first 
set of measures are reiterated, so to appreciate how the 
efficiency of the HSPN is recovered when inhabitants are 
reallocated. The results are presented in Figures [TT] [l2 
andfUl 



In order to compare the several reconstruction pro- 
cesses it is possible to consider the areas under the re- 
covery curve Y{C) = E{C) / E pre- shock- In engineering 
the integration of the area under the curve that gives 
the recovery of the quality of an infrastructure as been 



labelled resilience id 



R = 



t2-h ' 



where t is the time in post-event days and ti and t2 
are the endpoints of the time interval under consider- 
ation (Reed et al. 2007). Instead of counting the time 



steps after the event, we measured the progress in terms 
of reallocated people C. The introduction of a further 
normalized function 



yic) 



Y{C)-Y{0) 
1 - Y{0) ' 



which has value = 1 for = 1 and value ^ = 

in the aftermath Y(C) = 1^(0), allows to compare the 



^ Some authors refer to the resihence as the area over the recovery 
function ( jBruneau et al. [2003^. 
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FIG. 10: Identification of the new residential areas (in 
black circle) in the new sites recovery strategy 
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recovery processes which ensue after events of different 
magnitude. In formula: 



n 



J^-- y{C)dC 



Crr 



(9) 



where Cmax is the total number of people to be allocated 
after an earthquake of a given PGA. With this choice a 
resilience of 7^ = 1 corresponds to an instantaneous and 
total recover. 

The values of resilience 71 are reported in table [llj Ob- 
serving the results in terms of recovered efficiency for the 
different strategies, a number of remarks can be made. 
First of all it can be observed that the new sites strategies 
provide with slower efficiency recovery, if compared with 
the status quo strategies. In particular, in the new sites 
strategies after the first step, consisting in the realloca- 
tion of the 20% of the inhabitants into new buildings, the 
efficiency values present a negligible increase. Further- 
more, at the last step the efficiency values are not totally 
recovered to their initial values. This means that the new 
configuration of the city, with 4 new buildings used to re- 
allocate the 20% of the displaced citizens, is less efficient 
than its original configuration. A further remark can 
be made with regards to the up-down and down-up op- 
tions, corresponding to rehabilitating the buildings from 
the largest to the smallest and vice-versa, respectively. 
It can be observed that the most efficient option can not 
be unambiguously identified. 

One more remark can be made with regards to the 
status quo inwards and status quo outwards strategies. 
We can observed that for earthquake of PGA = 0.2 and 
PGA = 0.25 the status quo inwards option is the most 
resilient /fastest strategy in recovering the efficiency Ecc^ 
with respect to all the other strategies investigated. This 
result can be justified with the fact that, in the consid- 
ered case study, the more vulnerable buildings, i.e. the 
masonry buildings, are in the centre of the city. Hence, 
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FIG. 11: Efficiencies E^c for the 4 recovery strategies in 
the 0.2PGA(a), 0.25PGA(b) and 0.3PGA(c) 
earthquake scenarios. The absolute values are 
normalized with the pre-shock efficiency Ecc,pre- 
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TABLE II: Resiliencies IZ of the different reconstruction 
processes, computed using definition [9] For each couple 

HSPN and PGA, the rows highlighted in green 
correspond to the most resilient reconstruction strategy, 
while the rows highlighted in red correspond to the two 
less resihent strategies. 



(c) 

FIG. 12: Efficiencies Ecg for the 4 recovery strategies in 
the 0.2PGA(a), 0.2bPGA{h) and 0.3PGA(c) 
earthquake scenarios. The absolute values are 
normalized with the pre-shock efficiency Ecg^pre- 
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in the post event configuration the center of the city is 
much more disconnected, looking hke a "hole" in the 
HSPN networks (figure [9|. Thus, in the outwards strat- 
egy the first steps are not so efficient since the building 
are restored onto a seriously damaged network. On the 
contrary, in the inwards strategy, the restored buildings 
are progressively reinstalled onto a network which is also 
progressively reconnected, making this strategy more re- 
silient. 
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FIG. 13: Efficiencies Ecs for the 4 recovery strategies 
the 0.2PGA(a), 0.25PGA(b) and 0.3PGA(c) 
earthquake scenarios. The absolute values are 
normalized with the pre-shock efficiency Ecs,pre- 
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IV. CONCLUSIONS 

In this paper we proposed a novel methodology to 
quantify the resilience of complex social-physical urban 
systems against disasters. This methodology aims to 
merge together the two classical approaches to resilience: 
the engineering resilience, generally meant as the capa- 
bility of a system to recover its initial configuration after 
a shock, and the ecosystem resilience, generally meant 
as the capability of a system to recover its functionality 
after a shock, even in a new configuration of its compo- 
nents, different from the pre-shock one. The procedure 
presented here is inspired by the idea that city resilience 
should be assessed taking into account its social compo- 
nents, namely the citizens, which are the final users of the 
urban system as a whole. Our approach to quantify city 
resilience is based on the efficiency of hybrid networks 
composed by humans and urban infrastructures. In or- 
der to assess the capability of a city to recover its func- 
tionality after a shock event, we compare the efficiency of 
the corresponding hybrid networks before and after the 
shock event has occurred. We have discussed the results 
obtained on a case study, the city of Acerra, for which 
we simulated different earthquake scenarios and we anal- 
ysed three hybrid networks, namely the networks pro- 
viding residences, goods and schools. Different recovery 
strategies have been investigated, consisting in restoring 
the pre-event configuration or establishing a new config- 
uration, with different reconstruction priorities. 

While the main idea of the methodology presented here 
is to assess the resilience of urban recovery, the quantifi- 
cation of the efficiency of hybrid networks can be em- 
ployed also for other purposes, e.g.: 

1. To compare the efficiency of different urban con- 
figurations or different urban planning strategies. 
This would be a urban planning task; 

2. To design the reconstruction operations after an 
hazardous event; the best reconstruction strategy 
can be selected, identifying the physical configura- 
tion which maximises the performances of all the 
hybrid networks. 

It is underlined that, with regards to the assessment of 
post-event recovery strategies, the procedure here imple- 
mented does not account for the availability of financial 
resources. In fact, resilience is quantified in terms of re- 
allocated citizens and not in terms of elapsed time after 
the event. Thus, introducing the financial availability 
the time to reallocate citizens could be evaluated, that 
is the time to restore the damaged buildings and infras- 
tructures. However, the availability of financial resources 
can be considered a social-economic background input, 
not depending on the adopted recovery strategy. 

Once the financial "weight" is introduced, the proce- 
dure could be enriched considering different rehabilita- 
tion costs for different buildings and infrastructures. As 



an example, it is estimated that the cost of rehabilitating 
masonry buildings in the historical centre, in Italy, given 
the damage state, is twice as much the cost of rebuilding 
reinforced concrete structures; and construction time is 
also significantly higher. This would affect the perfor- 
mance of the outwards strategies that, starting rehabili- 
tating the historical centre, would be even less effective 
than the inwards strategies. 

Currently, further research activities on this topic are 
ongoing. The intention of the authors is to enrich and cal- 
ibrate the measures on hybrid networks, to achieve a set 
of quantifiable efficiency values, to be used as proxy mea- 
sures of the quality of life of the inhabitants of the urban 
systems. The final aim is to build a reliable tool to assess 
the quality of life in urban environments, in order to have 
an engineering quantification of the city resilience, with 
an ecosystemic and social-centric perspective. 
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